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Numerical studies of the flux creep in superconductors show that the distribution of the magnetic 
field at any stage of the creep process can be well described by the condition of spatial constancy of 
the activation energy U independently on the particular dependence of U on the field B and current 
j. This results from a self-organization of the creep process in the undercritical state j < j c related 
to a strong non-linearity of the flux motion. Using the spatial constancy of U, one can find the 
field profiles B(x), formulate a semi-analytical approach to the creep problem and generalize the 
logarithmic solution for flux creep, obtained for U = U(j), to the case of essential dependence of U 
on B. This approach is useful for the analysis of dynamic formation of an anomalous magnetization 
curve ("fishtail"). We analyze the quality of the logarithmic and generalized logarithmic approxima- 
tions and show that the latter predicts a maximum in the creep rate at short times, which has been 
observed experimentally. The vortex annihilation lines (or the sample edge for the case of remanent 
O ' state relaxation), where -8 = 0, cause instabilities (flux-flow regions) and modify or even destroy 

^' ( , the self-organization of flux creep in the whole sample. 
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I. INTRODUCTION 



Since the discovery of the giant vortex creep jl]] in high-temperature superconductors (HTSC), it has become clear 
that the relaxation processes in these compounds may be very rapid compared to usual low-temperature supercon- 
O | ductors. The magnetization current j and, in turn, the magnetic moment M, which is approximately proportional to 
j in most cases, drop considerably during the usual experimental time windows of a few hours (or even less) down to 
| small values j <C j c , in particular, at elevated temperatures. Here j c is the critical current, which divides the regimes 
^ . of flux creep (j < j c ) characterized by the Boltzmann factor exp(— U/kT), where U(B,j,T) is the activation energy 
' for flux creep, and the non-activational flux flow (j > j c ) with U = 0. Due to such a pronounced relaxation, both j 
and M are determined in HTSC mostly by the flux dynamics in contrast with conventional superconductors where 
relaxation is usually very slow, so j = j c for any accessible time windows. This has given rise to an extensive study, 



S-j I both theoretical and experimental, of magnetic relaxation and vortex dynamics in HTSC (see as reviews). Most 



of them are based on the logarithmic solution 0: 
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U~kTln[ 1 + 1 ) , (1) 



where to is the logarithmic time scale for flux creep. We will discuss Eq. (0) in detail in Section III. 

A closely related problem is the so-called "fishtail" effect, i.e., the anomalous increase of M as a function of H M, 
or the increase of locally measured j as a function of B |7j|§] . This effect serves as a test for different models of flux 
pinning and creep. The crucial question is whether the non-monotonous behavior of j results from the same feature 
in j c ("static fishtail"), or arises from a faster relaxation of j at small B, whereas j c is a monotonously decreasing 
function of B and itself shows no anomaly (" dynamic fishtail" ) . The latter possibility implies that the fishtail effect 
should disappear at shorter time windows or lower temperatures where the effect of relaxation is negligible and j = j c 



IB- 
After the instantaneous switching on of the external field H, the flux-flow process develops towards establishing 

a nearly critical profile j = j c {B), where the Lorentz force (j c X </>o) /c is compensated by the pinning force (<po is 

the flux quantum and c is the velocity of light) all over the sample. Usually the duration Tfi ow of flux flow does not 

exceed a few milliseconds (see Ref. |l0| and references therein). As j drops below j c , the slow process of flux creep 

starts. The creep rate is mostly determined by the Boltzmann factor exp(—U/kT), where U = U(B,j). Of course, 

U depends also on temperature, but the creep experiments are usually conducted at constant temperature. In many 

cases the dependence of U on B can be neglected, which implies a crucial simplification for the theoretical description 

of flux creep. For instance, in an infinite slab of width 2d (— d < x < d) in the parallel external field H, the variation 

of the magnetic induction SB = B(d) — B(0) do not exceed H* , where H* = (47r/c)j c d is the full penetration field IO] 
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(j c is considered to be field-independent). For H ^> H* one can neglect SB < H* . If both j c and d are sufficiently 
small, say, j c ~ 10 4 A/cm and d ~ 0.1 mm, then if* is of order hundreds of Gauss. In this case the above condition 
is easily fulfilled for most if, and the activation energy appears to be field-independent: U(j, B) = U (j, if). Then the 
field profiles B{x) are almost straight |l^,|l^], i.e., j = const throughout the sample at all creep stages. 

However, in larger samples (d> 1 mm) with strong pinning {j c > I0 5 A/cm 2 ), one gets SB ~ if* > 1 Tesla, which 
implies that the spatial variation of B may not be small compared to H . Of course this estimation may not be valid 
at high temperatures where j c drops. But for large and not too clean samples well below T c the dependence of U on 
B is essential, and the field profiles are not straight. On the other hand, due to very fast relaxation B can vary by 
orders of magnitude during experimental time windows of order of hours (see, for instance, ffl]), which also requires 
taking the [/(independence into account for the consistent description of the relaxation process, especially in the 
dynamic models of the fishtail formation. 

The goal of this paper is to study the relaxation process for various dependencies of U on j and B. In Section 
II we show that the field profile at any stage of the relaxation process can be described by the condition of spatial 
constancy of the activation energy: 

U(x) ^ const (2) 

throughout the whole sample, where const depends on time only. Since U — U(B,j), Eq. (Q) provides an implicit 
relationship between B and j, manifesting a condition of self-organization of flux motion in the undercritical regime 
j < j c . In other words, according to Eq. (||) the field profiles form a one-parameter family Bu{x). The problem to be 
solved in order to describe the flux creep is to find these profiles together with the dependence of U on time. Note 
that this case differs from the self-organized criticality (see the pioneering work |l3j and further applications to the 
superconductors in a critical state |14[| ), where the peculiarities of the flux motion (avalanches, critical exponents, 
etc.) are considered in the vicinity of the critical state j = j c (he., at U = 0). The condition j = const found in 
previous studies is obviously just a particular case of Eq. (||) provided U is independent of B. In Section III we 

analyze numerically and, using Eq. (^]), semi-analytically the flux creep for various U(B, ^-dependencies, particularly 
for the most general collective creep behavior U oc B a j~^. We show that at short, but experimentally available time 
scales the creep process differs significantly from the logarithmic solution (see Eq. ([!])) and shows a maximum in the 
relaxation rate, dU/dlnt, in accordance with the experimental data. The semi-analytical solution provides a good 
approximation to the exact (numerical) one at all time scales. In Section IV we apply these ideas to the problem of 
the anomalous magnetization curve (fishtail) and show how the dynamic development of the anomaly in j (or the 
same, in M oc j) can be described semi-analytically. In Section V we study the effect of so-called "annihilation lines" 
in infinitely long samples, where B changes sign and vortices with opposite directions annihilate each other, on the 
self-organization of the flux motion. A particular case of such a line is the edge of the infinitely long sample in the 
remanent state. The vortex velocity v shows a peculiarity (divergence) at such a line, resulting in the appearance 
of flux-flow regions in the vicinity of the annihilation lines. We show that these peculiarities affect deeply the self- 
organization of creep and the condition of spatial constancy for U (see Eq. (g)) is modified or even destroyed in the 
whole sample (and not only in the vicinity of the annihilation lines). 



II. SPATIAL CONSTANCY OF U 



Consider an infinite slab of thickness 2d (— d < x < d) with the magnetic field B(x) parallel to z-direction and the 
current j flowing along y. The external field H is switched on instantaneously when no vortices are present in the 
sample, which corresponds to zero-field-cooled experiments. We consider H ^> H c i, where H c i is the lower critical 
field, and disregard the effects related to the latter. 

The flux creep is described by the diffusion equation Jl5|] : 

dB JD 

dt ~ dx' [6) 

where 

D = Bv = A—Bj exp(-U/kT) (4) 
cn 

is the magnetic flux current, v is the vortex velocity, n is the Bardeen-Stephen drag (friction) coefficient |l6| ] for flux 

flow and A is a numerical factor. Note that D is proportional to the electric field E — (B x v)/c in the sample. The 
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form of the magnetic flux current D is chosen such that at U — and A — 1 the flux velocity v corresponds to the 
Bardeen-Stephen expression fl6|| for the flux flow: v — vn ow = (j>oj/cn. 

It has been already discussed [E3 that the strongly non-linear Eq. (^) should obey a self-organized behavior. This 
means that if a fluctuation SU appears in the sample, it results in a significant (exponential) local change of the flux 
current D cx exp(—5U/kT) which, in turn, leads to fast smearing out of the fluctuation. In other words, \SU\ ~ kT 
is the scale of "permitted" variations of U. This has been proved experimentally by direct measurements of U using 
the Hall probe technique jf7|. We suggest a more general criterion of self-organization as follows. 

The variation of the flux current density D (or the same, of the electric field E) within the sample can be written 
as: 

where B, j and U are formally considered as three relaxing parameters, though U = U(B,j). If one of the three 
terms in the right-hand-side of Eq. (|^) appears to be considerably greater by absolute value than the other two terms, 
the corresponding parameter governs the relaxation process, i.e., relaxes towards its mean value irrespective of what 
happens with the other two. Obviously, this leads to a self-organization of the flux diffusion process which implies 
that the three terms in the right-hand-side of Eq. (^) tend to keep the same order: 

3D 3D dD ... 
dB 5B = ^J 5] = dU 6U - (6) 

The above condition can be considered as a mutual confinement for variations of B, j and U. Taking the expression 
for D from Eq. (||), we get a limitation for SU: 

SU (SB Sj\ 

— <max ■ ( 7 ) 



kT ~ ' { B j 

Note that the above estimation does not require the condition U kT, i.e., it should hold starting from the very 
early stages of flux creep. 

If one of the three terms in Eq. (^J) is very small (or absent) for any "external" reason, then the self-organization 
applies to the two other ones. For instance, if B is much greater than H* and thus SB/B -C 1, we get: 

SU ^ dUSj_ ^ UoS_l 

kT dj kT kT j c U 

where U c is the characteristic activation energy for j — > 0. Since in general U c S> kT, we get Sj < (kT/U c ) j c <C j c , 
as has already been discussed in Ref . PJ . 

It is worth mentioning, however, that at the locations where j = or B = the variations of U can exceed kT 
significantly, as follows from Eq. (^) . The first of these two conditions (J = 0) regularly holds at the center of the 
sample, and we will comment on this point in the following Section III. The second condition (B = 0) holds at the 
lines where vortices of different sign annihilate each other, or just at the edge of the sample in the remanent state. 
We will devote special Section V for the latter case. 

Eqs. (Q)-(@) prove the spatial constancy of U throughout the sample with a kT precision (see Eq. @). The analytical 
results based on Eq. (0) we will refer below as "semi-analytical" ones. 



III. ONE-DIMENSIONAL CREEP EQUATION 

In an infinite slab the current j is related to B by the Maxwell law: 

c dB 

> = -^ (9) 

if one uses a reference system xyz where B || z and j \\ y. For a platelet sample in a perpendicular external field, 
where the in-plane field component B x appears, the relation between j and B becomes more complicated ^|,[l8| . Here 
we focus on the one-dimensional creep problem where B x — 0, so B z = B. 

After substituting Eq. (|^) into Eq. (g) one gets the basic one-dimensional equation for flux motion: 
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The numerical coefficient A depends on the creep mechanism and should not necessarily be of order of unity |Tg| ] 
(see also discussion in Rcf. |fio|j ). However, the low field measurements (i.e., in the single vortex pinning regime) in 
YBa2Cu307_ x crystals ]l7| , where Eq. ([!(]) was experimentally studied by direct local measurements of B, dB/dx 
and dB/dt, revealed that A Thus, at least in the case of single vortex pinning, Eq. ( |i"o|) is consistent with the 
equation for flux flow: 



dB _ d 

dt dx \4tt7] dx 



o B 8B 



(11) 



since the latter can be obtained from Eq. ( |To| ) at U — and A= 1. 

The features of self-organized criticality in the solution of Eq. ([h]) have been analyzed (20| for the case of switching 
on of a small additional field SB on the background of B ^S> SB already present in the sample, and for a specific 
(logarithmic) dependence of the activation energy on the current: U = Uq ln(j'o/j). In terms of the energy distribution 
U{x) across the sample the case considered in Ref. |20| implies that in the beginning the energy was very large (or 
infinite) in the whole sample, since j = 0, then an area of small U appeared at the edge (after switching on SB), and 
the propagation of this "fluctuation" of the U(x) profile was studied. 

In contrast with Ref. pof we consider below the instantaneous switching on or removal of the whole external field 
H. A sort of self-organization, i.e., establishing of a "partial critical state" [|TJ with j(B) oc j c (B) has already been 
reported for this case. We show below that this result follows from our more general approach based on Eq. (|^) 
if U = U(j /j c (B)), but for an arbitrary U(B, j)-dependence the partial critical state may not be established. Our 
general results on self-organization of the flux creep do not depend on the specific U(B, j)-dependence, but focus on 
the collective creep behavior U(B,j) oc B a j~ fl as mostly relevant in HTSC. We will not consider the time-dependent 
boundary conditions H = H (t). Some results for the latter case can be found in Refs. p^J^ . 

The integration of Eq. (10) can be performed as follows. Defining the magnetization as: 



1 

2d 



(B - H) dx 



and integrating Eq. ( |l0| ) over x, we get: 



dm _ fa 
dt Airnd 



H 



dB 



dx 



aq>(-U edg e/kT), 



(12) 



(13) 



X — ±rf 



where U e d ge — U(x — ±d) is the activation energy at the edges of the slab. For the case (discussed below) of straight 
field profiles: \dB/dx\ = const, and constancy of the activation energy (the latter implies that U e d ge can be substituted 
by the mean and almost constant U over the sample, see Eq. (||)), one can rewrite Eq. (JlJ) in the form: 



dm . m , ,, 

^ A- exp (~U kT , 
dt t 



(14) 



where t = 2irr)d 2 / (t>$H . 



A. Flux flow (j > j c ) 

After switching on the external field H , the flux flow starts and lasts until the vortices fill the sample up to the 
critical profile j c . Its duration Tfi ow can be easily estimated if one notices that already during the flow regime 
the B(a;)-profiles are almost straight, i.e., \j\ = (c/Att) \dB/dx\ ~ const (see dashed lines in Fig. la). Using the 
straightness of the field profiles, one can estimate the magnetization m as: 

for j>j*, (15) 

for j<f, (16) 

where the current j* — cH/Aird shown as a grey solid line in Fig. la discriminates between the incomplete and 
complete penetration of flux into the sample. Note that we have chosen H > H* , where H* — 4irj c d/c is the field 
of full penetration. This means that j c < j*, i.e., by the completion of the flux flow stage, flux penetrates the whole 
sample (see Fig. la). Then, substituting these values into Eq. ( |l3| ) and solving it at U = and A = 1, we get: 



m{j) = H- 



cH 2 
8njd 



m(j) ^ 



2njd 
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j = f]/Y t > t<T/2 (17) 

jSj-expQ-iV t>r/2 (18) 

where r/2 = mjd 2 /4>qH is the time of full penetration in the flux flow regime. For j c < j* the time of establishing 
of the critical profile after switching on the external field is Tfi ow = t(j c ) — (r/2)(l + 2 In j*/j c ) — t as follows from 
Eq. @. 

The crossover from the flux flow to flux creep process is well defined, i.e., the critical profile j — j c is established 
at t = Tj\ ow almost exactly throughout the whole sample (see Fig. la). More exactly, the fluctuations of U which 
appear in the whole slab at the crossover from flow to creep (j = j c ) are of order kT according to Eq. (|^), therefore 
Sj < (kT/U c )j c <C j c - The only exclusion is the very center of the slab, x = 0, where j = and U shows relatively 
strong variations, as we discussed in the previous Section. However, this area is narrow and can be neglected when 
considering the profiles B(x) and magnetization m. 

After the flux flow stage is completed, a much slower process of flux creep starts, and various cases of U(B,j) can 
be analyzed. First we consider the simplest case where U depends only on j. 



B. Creep at j c — const, U = U(j) 

This case has been already studied in Refs. [^0| , ^2|j2l]j2^| . Here we analyze it as a test for our numerical solution 
before consideration of more complicate models of U(B,j). During the stage of flux creep (j < j c ) (see Fig. la), the 
field profiles are even more straight than during the flux flow stage [ ^0|Jl2| , i.e., \ j\ = const, and, in turn, U(j) = const 
(see Fig. lb). Note a very narrow increase of 5U at x = 0, where j = 0, which is consistent with the comment at the 
end of Section II. 

Since U e dge is approximately equal to the mean U over the sample, one gets from Eqs. (|l3|) and ([H^), using also 
Eq. (§): 



dll _ dU dj dm ^ A . 
dt dj dm dt r 



dU 



dj 



exp(-U/kT). (19) 



This equation can be integrated numerically for any [/(j)-dependence, and also can be solved with a logarithmic 
accuracy ||, see Eq. (|]). The latter means that the real [/(j)-dependence is substituted by the tangent straight line 
with a slope dll/dj, as shown in Fig. 2, which is reasonable since the relaxation slows down exponentially as U grows. 
Thus the system spends most of the relaxation time near the final point where U(j) and its tangent line almost 
coincide. Such an approximate solution of Eq. (|^|) acquires the form of Eq. (^) with 

kT 2-nnd 2 _ rkT 
t0 ~ \dU/dj\A<f>ojH ~ Aj\dU/dj\- ( } 

As becomes clear from Eq. (lj) and Fig. 2, to is the time required to get from U = — oo (which corresponds to t = —to) 
to U = (which corresponds to t — 0) along the non-physical part of the tangent line, corresponding to negative U . 
Thus to has no direct physical meaning and should not be mixed with the characteristic duration Tfi ow ~ r of flux 
flow, see Eqs. (R-©. 

Eqs. (|l|) and (|20[) provide a logarithmic approximation for the time required for a system to reach the energy U. 
However, in order to use Eq. ([j]) to describe the [/(t)-dependence one observes that to is not actually a constant and 
depends on U (or the same, on t). This effect is not of great importance at dU j dj ~ const, i.e., where U is an almost 
linear function of j, but cannot be neglected in the opposite case of strongly non-linear dependence of U on j, where 
dU /dj changes significantly. 

Consider the case of the collective creep model 

U = U C [{j c /jY - 1] , (21) 

which is an example of such a non- linear dependence. Here the exponent fj, varies Q from /i = 1/7 (single vortex 
creep) to /i — 5/2 (small bundles). The straightforward solution of Eq. (Il9]) then gives: 



5 



where Ei is the integral exponential function. The logarithmic approximation for Eq. ( |22| ) acquires the form: 

U = kThx(l + AfJ,-(U + U c )/kT] , (23) 



which, if compared with Eq. (Q), implies: 

TkT 



Afi(U + U c ) ' 



(24) 



If during the creep process j decreases down to j mm <C j c , then the energy increases up to U max ^> U c , where 
[/ max = [/(jmin)j an d to decreases down to t^ n = rkT j AfJ.U max . Regularly, the latter estimation, to = t™ m , is 
substituted into the logarithmic solution Eq. ([j]), and time dependence of to is neglected. This results in an almost 
linear dependence of U on ln(i) (see Fig. 3). Let us call such an approximation (where to is treated as a constant) as 
a "pure" logarithmic solution, whereas Eqs. (p3|)-(p4]) provide a "generalized" one. 

The straightforward solution and both generalized and pure logarithmic ones are compared in Fig. 3. In the same 
figure we show an exact numerical solution of Eq. ( flO| ) obtained without any assumptions on spatial constancy of U. 
We used r / A[i as a useful time constant in Fig. 3. One observes that the generalized logarithmic solution works at all 
t, and together with the straightforward solution provide a perfect fit to the exact one. On the other side, the pure 
logarithmic solution shows significant deviations from the exact one, especially at short times. This is consistent with 
Eq. (p4|), since there U < U c , and, in turn, to ^ ^o" n - Particularly, the logarithmic solution misses the characteristic 
maximum in the creep rate dU/dhit which appears at Afit/r = 10. For n ~ 10~ 5 g/cm-sec (which implies that 
T is well below T c ), A d 1, fi = 1/7 (single vortex creep) and d = 1mm, the maximum in dU/dhit appears at 
t* ~ 60 /H sec, where H is measured in Oe. Thus the maximum or at least its tail can be resolved at experimentally 
accessible times. Note that the numerical factor in the latter estimation grows oc d 2 , and for some larger samples t* can 
be significantly greater. The position of such a maximum determined experimentally can be used for determination 
of the parameters t/A and /1. 

The above theoretical results are compared in Fig. 4 with the experimental data obtained on a large (d = 1 mm) 
YBa2Cu307 crystal. The relaxation of the magnetic moment was fitted by the direct numerical solution of Eqs. (|l4| ) 
and (pi]), where t/A, j c , H and U c were considered as independent fitting parameters. Thus [/(t)-dependence is 
found directly from the experiment and is compared in Fig. 4 with the theoretical curve found from Eq. ( p2| ) with 
the same parameters. The experimental and theoretical results almost coincide, and the characteristic maximum in 
d(U/kT)/d\nt at t = lOr/An is very clear. 

C. Creep at j c = const, U = U {B/B ) (j c /j - 1) 

Consider this simplest model dependence of U on B and j, which can mostly be analyzed analytically before a more 
complicated case of collective creep. We have changed the notation Uo instead of U c , as in the previous subsections, 
for reasons which are clarified below. Fig. 5 illustrates the numerical solution for B{x) and U(x) profiles in the case. 
One observes that the general condition U = const holds in this case as well as in the previous one, where U was 
independent of B. The same narrow peak in U is located at the center of the sample x = 0, as described in Section II. 
Since U depends only on the ratio Vol Bo, we can choose Bo arbitrarily. It is most convenient to accept Bq = H*. 
Then, using the Maxwell equation (^) and the condition U = const, one obtains the approximate expression for the 
field profile in the sample: 

- = 1H 1 In—. (25) 

d H* Uo H v ' 



Denoting the field at the center of the sample as B(0), which is found from Eq. (25) at x = 0, one can rewrite the 
magnetic moment m (see Eq. (|l2|)), as: 



rn 



H 

xdB. (26) 

B(0) 



Then 



dm dm dU 1 
~dt = dU~dt = ~Z/o 



dU 
~dt' 



G 



On the other hand, as follows from Eq. (|l3|), 



dm A 
~dt = ~4r 



dB 



dx 



exp(-U/kT), (27) 

x—znd 



and one gets an equation which determines the activation energy: 

dU AU h 



dt t 2(h + u)[h-b(0)-b(0)ln(h/b(0))] 



exp(-Z7/fcT), (28) 



where we denoted h = H/H*, b = B/H*, and u — U/Uq. This cumbersome expression can be reduced if one uses the 
expansion over A = [H — B(0)] /H, which becomes exact at A — » 0, but actually holds with 10% accuracy at worst 
for all A < 1, i.e., for all B(0) < H. Then one gets: 

^L^A {U + Uc) (i + ^a) eM-U/kT), (29) 



dt t v ' \ 3 

where U c = hUo. This result differs from Eq. ( p3| ) only by the correction factor 1 + 2A/3. This means that the results 
of the previous subsection, where U was field independent, apply here just by accounting for the correction factor 
1 + 2A/3, which in most cases is not of great importance. 

Eq. ( p9| ) enables us to find the [/(i)-dependence, and then, using Eq. ( |25| ) we get the field profile B(x) as a function 
of time, i.e., solve the creep problem completely. We call such an approach as "semi-analytical" solution. Its results are 
compared in Fig. 5a with the exact solution obtained by a direct numerical integration of Eq. (|l^) with no assumptions 
on constancy of U. One observed that the semi-analytical solution, being much less time-consuming (note that the 
solution of Eq. ( p9| ) is quite universal and has been already obtained in the previous subsection) , provides a perfect 
fit to the exact description of the creep process. 



D. Creep at j c = const, U = Uo(B / Bo) a {jc/j) 1 * (collective creep) 

This is the general dependence of U on B and j in the collective creep theory [|| for j -C j c - The spatial constancy 
of U holds in this case as well as in previously considered cases. We have checked it numerically for various a and fi. 
The condition U = const together with the Maxwell equation (|p) determines the field profile: 



W 



— ( 



1 



(30) 



where we denote v = 1 — a/fx and assume Bq — H* , as in the previous case. Here we take a/(i (the case a — Li is 
almost identical to that considered in the previous subsection). The field 6(0) is determined, as follows from Eq. (|3C|), 
by h" — b(0) v — v/u 1 ' 11 . The magnetic moment can be calculated using Eq. (126 



H* 



and, using Eq. (pTj), one gets: 

dU _ Au{v + l) 
~dt ~ 2rh' / - 1 u 2 /v 



6(0) (h v + u- 1/f *)} Uexp{-U/kT) 



(31) 



(32) 



This expression, being a bit cumbersome, can be reduced using the expansion over A = (H — B(0))/H, which, as in 
the previous case, works with reasonable accuracy (better than 10%) at all B(Q) < H: 



dU 
~dt 



Au 



U 



1 — A 

3/i 



exp(-U/kT). 



(33) 



This result is very similar to Eq. (|29j). The absence of U c in Eq. ( |33| ) corresponds to the absence of term —1 in the 
dependence of the activation energy U on j in this model (compared with the two previous cases). The factor in the 
square brackets in Eq. (33) describes the effective renormalization of U in the pre-exponential factor resulting from 
the dependence of U on B. If a = 0, which means that U is independent of B, then the renormalization disappears. 
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The same happens at B(0) —* H, At a = \i the correction factor reduces to 1 + 2A/3, which is consistent with the 
previous case, where a = (X = 1. 

In Fig. 6 we compare the direct numerical solution of Eq. ( |To| ) with the semi-analytical one determined by Eqs. (|3(i|)- 
(|3l|). Fig. 6a shows the numerical (exact) B(x) profiles compared with Eq. (^o|), and Fig. 6b shows m(lni), obtained 
numerically from Eq. ( |l0| ) and semi-analytically from Eqs. (|o|)-(^3|). The quality of the semi- analytical approach is 
perfect in this case as well as in the previous one. 

In the most general case: 

U = U Q (B/B a ) a [{j c /jf - 1] (34) 
one gets an expression which naturally conforms to Eqs. (^9|) and (|33|): 



at t 



3u 



exp(-U/kT), (35) 



where, as above, U c = KUq. We skip the cumbersome derivation of the last expression, which requires expansion over 
A starting from the equation for the field profile B(x). The generalized logarithmic solution of this equation acquires 
the form 



U = kTln ( I + + U °] P [ + WW^t ) , (36) 

2irrid /! kT 



which coincides with Eq. (Q) if 



tQ ~ AfoH^U + U c ) [1 + (2a/3 M )A] ~ IT H ' ( ' 



where we introduced tn = to(H — H*). Note that tn is almost field independent, since H enters tn only via A. Note 
that both to and tn depend on time via U and the correction term (in square brackets) . 



E. Creep at j c = jc(B) ^ const 

Above, we have considered only j c = const. However, the field dependence of the critical current, j c = j c (B), does 
not violate the general condition, Eq. (|^), of self-organization of flux creep. As a direct consequence of this condition 
it is worth mentioning the following: If the dependence of U on j and B has the form: U cx f(J c (B)/j), where / 
is an arbitrary function, then the spatial constancy of U results in establishing of a "partial" critical state (2lj with 
j cx jc(B). For instance, if the critical current obeys the Kim dependence j c (B) = joBo/(Bo + B), then the field 
profile B(x) during the relaxation should be determined by the condition j = pjoBo/ (Bq + B) with < p < 1. 

However, for more complicated dependencies of U on B and j this is not the case, and the profiles B(x) can differ 
significantly from that in the critical state. In the next Section we consider an example of such a behavior. 



IV. SEMI- ANALYTICAL SOLUTIONS FOR ANOMALOUS MAGNETIZATION (FISHTAIL). 

The equation ( |35| ) and its reduced forms (see Eqs. (^9|) and fl33])), which are just ordinary differential equations, 
present the method of semi-analytical integration of the equation for flux motion (see Eq. (|Io|)), for the case of 
collective creep, where the dependence of U on B and j is described by Eq. (0), or by its reduced versions. Of 
course an analogous solution can be found for any U(B, j)-dependence, not only for that described by Eq. ( pih. T his 
semi-analytical approach provides a good fit to the exact solution, obtained by numerical integration of Eq. (|10|), as 
one can see in Fig. 6. The correction factor 1 + (2a/3/x)A can be neglected except for short times t > t/Au. 

The semi- analytical solutions can be applied for the description of an anomalous magnetization, coined a "fishtail" , 
found in clean high-T c superconductors |i|-f§. Note that j c enters Eq. ([55]) only via the correction factor which is 
negligible in most cases, especially at high fields H H* where A « 1. Thus the solution U(t) of Eq. ( |35| ) is 
determined by the current exponent /i (and not by the field one, a), by the characteristic energy U c , and by r (which 
in turn depends on d, r\ and H). If one measures the magnetization current j at the edge of the sample, where 
U(B,j) = U(H,j) as a function of H, keeping the time window t of the experiment constant for each H (this is the 
case for most studies of fishtails), then U along the measured line M(H) or j{H) can be written, as follows from 
Eq. (||), as 
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U(H, t),kT - H U + U c ) ^lnH + tat + In (^^^) , (38) 

where the last term in Eq. (|3^ ) almost docs not depend on H and t. This means that the magnetization curve is 
determined by: 

dU dU _ dH 

kT~ u + u c ~ h ' { ' 

In Fig. 7 we present the results of our semi-analytical approach to the problem of dynamic fishtail formation taking 
j c (B) = joB c /(B c + B) (Kim model) and collective creep with 



U = U {B/H* 



MB) 



J 



- 1 



B c + B 
B 



(40) 



The last factor in this equation is added to cancel the dependence of j c on B and keep the general collective creep 
condition: U oc B a j~^ for j -C j c - At each H we find the energy U down to where the system relaxes during the 
"experimental" time window t, and then, using this U, we determine the corresponding j x =±d according to Eq. fl4C|). 
The results show a clear fishtail due to fast relaxation at low fields (see Fig. 7). 

Note that Eq. fllpf ) provides an example of the case where the field profiles B(x) are significantly different from the 
critical one at j = j c , and a "partial critical state" ^ij is not established. 

Since dU = (8U/dH)dH + (dU / dj)dj , one obtains using Eq. ( |39| ) that the magnetization curve (fishtail) is deter- 
mined by the condition: 



dj [dU kT U + U c \ [dU 



dH \dH H U + U C ~ kT J \dj 



(41) 



For the case of collective creep, where U oc H a j M , we have dU/dH = (a/H)U and dU /dj = —(fj,/j)U. Then, taking 
into account that U + U C ^> kT, we get 

*'*?£fl-~V (42) 



dH fiH \ a U 

The peak of the fishtail, where j(H) reaches maximum, corresponds to U = kT/a, as follows from Eq. (fl2|). This 
implies that j increases as a function of H until it almost reaches the j c (H) curve. Far below j c , where kT/aU <C 1, 
one gets from Eq. Q t hat j oc H a ^. 

We see from Eqs. (|38|) and (^) that U changes along the magnetization curve obtained at a fixed time window 
t. However, one can measure jnt(H) keeping the product Ht as constant, which, according to Eq. ((3^), should 
result in a constant U along the magnetization curve (neglecting the correction factor 1 + (2a/3/i)A). The difference 
jm{H) — j(H), where j(H) is taken at t = const, is determined by: 

d{jm - j) kT 

dH H(dU/Bj)' [ ' 

which provides a tool for independent analysis of [/(j)-curve. 

Above in this Section we have considered the exponents a and n to be constants. However, different regions in the 
j — H diagram correspond to different relaxation regimes, such as single vortex creep, small and large bundle creep, 
etc. (see 0-0]). The energy scale Uo, as well as the exponents a and fj, may vary significantly from one region of 
j — H to another. As one observes from Eq. (^6|), the crucial exponent of the above two is \i. Its rapid change at the 
boundary between the creep regions from jii to /12 is equivalent to a change of H by factor ■ As follows from 

Eq. (|3q), this results in a change of ~ ln^a/Mi) m U/kT at the boundary between two creep regions. Thus U does 
not change much at the crossover from one pinning regime to another. However, j (and, in turn, M) can be changed 
significantly at such a boundary, since for different relaxation laws (different U c , a and fj,) the same U is reached at 
significantly different j. As H increases, the growing vortex bundles lead to increase of characteristic energies U c , 
thus one should expect a step-like increase of j when crossing the boundaries single vortex pinning — > small bundles 
— * large bundles. 

If one measures the exponent fx along the magnetization curve (see, for instance, Ref. 0), then a curve of constant 
U in the H — j diagram can be plotted using rather t oc (^tif) -1 instead of t oc H , as was suggested above. 
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V. RELAXATION IN THE REMANENT STATE AND ANNIHILATION LINES 

A particular and very interesting case, where the discussed above self-organization of the flux motion should be 
modified significantly, is relaxation in the presence of annihilation lines B — 0. The vortices and antivortices approach 
the annihilation line from different sides and annihilate each other. The arguments of Section II for the constancy of 
U are not valid in this case, at least in the vicinity of the annihilation lines (see comment at the end of Section II). 
Therefore, this case should be studied separately. 

Consider the simplest situation of remanent relaxation, where the field has been ramped up and then instantaneously 
removed, so B = at the edges x — ±d of the slab. There are no antivortices in this case, since the annihilation line 
coincides with the edge of the sample. 

The description of the flux motion in this case using Eq. (||) looks self-contradictory since at the sample edge 
B = 0, whereas the magnetic flux current D = Bv is finite at the edges and, moreover, obviously should reach there 
its maximal value over the sample. However, the contradiction is void provided the field vanishes at the sample edge 
as B (x y/d — x, i.e., proportional to the square root of the distance to the edge (see Fig. 8). At the same time the 
vortex velocity diverges at the sample edge as v oc dB jdx oc 1/ \Jd — x. This divergency is removed by an appropriate 
cut-off for d — x, which we discuss later in this Section, but inevitably leads to the appearance of the flux-flow region 
near the edge or, most generally, near the annihilation line. However D cx B{dB / dx) remains finite and continuous 
with no singularity at the edge. This is confirmed by direct computer simulations of the relaxation in the remanent 
state (see Fig. 8). 

Let us estimate the coefficient k in the square-root dependence B e d ge — k\J d — x near the sample edge, which 
should include the flux flow region {7 = 0. The magnetic flux current reaches at x = d its maximum over the sample: 
D x= d = (fa/inri) [B dB / dx] x ^ d — <j>ok 2 / 87777, but remains of the same order as the mean flux current (D) over the 
sample: D x= d — C (D), where C > 1 is a numerical factor. Estimating (D) ~ (Afa / Anrf) (B) (dB/dx) exp(— (U) /kT), 
(B) ~ B(0)/2, (dB/dx) ~ B(0)/d, one gets k ~ B(0)y/CAcxp(- (U) /kT)/d. Note that the above estimation is 
based on the constancy of U, i.e., U = (U) throughout the sample except the edge flux-flow regions, which we assume 
to be small. Thus we get: 



B edge (x) = 5(0)^04 cxp(- (U) lkT) d -^- (44) 

A natural cut-off for the area of applicability of Eq. (E3) is: d — x > A, otherwise the surface effects such as Bean- 
Livingston interaction with the surface p6| should be accounted for. There are additional restrictions: (i) the current 
cannot exceed the depairing one: j = (c/4ir)dB/dx < jd\ and (fi) the intervortex distance a ~ fa/ B(x) should not 
exceed the distance to the surface d — x at any x. It can be easily confirmed that the condition d — x > A is stronger 
than the other two at most reasonable values of B(0) and d. 

Let us call the region x < \x\ < d near the sample edges, where the activation energy grows from U = at the 
very edge (flux-flow region) up to U(x) = (U), as the area of "annihilation dominated" organization of flux creep. 
Its width d — x can be estimated as follows: We substitute Eq. (Q) into the collective creep formula for U(B,j), see 
Eq. (HU), and find x where U reaches its mean value (£/). This is of course a crude approximation, since Eq. ( f44| ) is 
valid in the flux-flow region only, and for the whole "annihilation dominated" region it provides an underestimation 
for B and, in turn, j. After straightforward calculations we get: 

d-x f(U)\^ f a-n(U) \ 
d -\U ) eXp [a + fxkT) 

The above result implies that the width d — x of the " annihilation dominated" region is crucially dependent upon 
the relationship between the exponents a and /i. For fi > a and (U) 3> kT this region appears to be exponentially 
small, i.e., x = d. Computer simulations show a step-like increase of U at the edge to the value comparable with 
(U), and then U grows smoothly and slowly towards the center of the sample (see Fig. 9a). Though SU appears 
to be significantly greater than for the case of finite H, discussed in previous Sections, even here U does not vary 
significantly: SU < AkT in the whole sample, excluding the sharp step at the edge. For the opposite case, fi < a, 
one finds from Eq. ( ff5| ) the unphysical result that (d— x)/d is exponentially large though, of course, (d — x)/d < 1 
anyway. This implies that our assumption about the spatial constancy of U(x) = (U) throughout almost the whole 
sample (except small edge regions) is self-contradictory in this case. Thus for /1 < a the effect of the annihilation 
line spreads over the whole sample, and there is not any evidence of constancy of U. This is confirmed by numerical 
simulations (see Fig. 9b). The boundary case, where fx = a, is illustrated in Fig. 9c. 
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(45) 



VI. CONCLUSION 



We considered the generalization of the logarithmic solution Eq. ([[]) for flux creep at different dependencies of 
the activation energy U on field B and current j and confirmed it by numerical analysis. The general condition 
which governs the relaxation is U(x) = const throughout the sample, and this result holds at any particular U(B,j)- 
dependence. This results from a self-organization of flux creep in the undercritical state j < j c , which implies that the 
influence of all the creep parameters, B, j and U, on the relaxation rate should be of the same order of magnitude. 
This self-organization should not be mixed with the self-organized criticality of flux motion at j = j c . 

For U independent of B, i.e., U = U (j), we restore the known result of straightness (j ~ const) of the field profiles 
throughout the sample. For the case where U essentially depends on B the condition of spatial constancy of U(B,j) 
determines the one-parameter family of field profiles Bjj(x) and enables us to find a semi-analytical solution for 
U(t) and, in turn, for time evolution of the field profiles Bu(x), i.e., to solve the creep problem completely. Such 
a semi-analytical solution provides a perfect fit to the exact numerical solution (obtained without any assumptions 
on constancy of U) and appears to be quite useful for the description of the dynamic development of anomalous 
magnetization (fishtail) due to fast relaxation rates at low fields. 

The effect of the annihilation lines B = on the self-organization of the collective creep where U cx B a j~^ (see 
Eq. (pi])) is crucially dependent on the relationship between a and fi. At a < u the effect is just an increase of 
variation 8U over the sample, with a step-like vanishing of U in the very narrow regions of flux flow in the vicinity 
of an annihilation line. However, for U ^> kT we still get SU -C U, i.e., the condition U = const holds qualitatively 
in this case. At a > /x the above condition no longer holds, and the presence of an annihilation line destroys the 
self-organization in the whole sample irrespective of it size. 
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Figure captions 

Fig. 1 a) The field profiles B(x) for flux flow (dashed lines) and creep (squares) for U = Uq [(j c /j) — 1] with H* = H/2. 
The critical state (j = j c ) and the full penetration state (j — j*) are shown in solid black and grey lines, 
respectively. Note an almost exact formation of a critical state at the crossover flow^creep. 
b) Spatial dependence of U at different times t/rfi ow . Note almost spatial constancy of U except narrow regions 
x = 0. 

Fig. 2 Relaxation of j can be imagined as a motion of a dot along the [/(j)-curve (solid). For the logarithmic approxi- 
mation U(j) is substituted by its tangent line (dashed). Time to corresponds to the "motion" along the negative 
part (U < 0) of the line, whereas Tfi ow corresponds to the "motion" from j = oo down to j — j c along U = 0. 

Fig. 3 Comparison of " pure" logarithmic (triangles) , straightforward (circles) and generalized logarithmic (solid line) 
solutions of Eq. ([L9|), which were derived under assumption U(x) = const, together with the "exact" numerical 
solution (squares) of Eq. ( |i"o| ) for U(j) = U c [{j c /jY — !]• Open and filled symbols correspond to U/kT and 
d(U/kT)/d\nt, respectively. Note that all the solutions except the pure logarithmic one show a maximum in 
dU/ dixit at Afit/r > 1. 

Fig. 4 Experimentally obtained relaxation rate d(U / kT) j d\nt (triangles) and normalized magnetic moment M/M c , 
which is equal to j / j c (circles) vs. t and their fit by the generalized logarithmic solution at the same values of 
parameters: A/t = 0.03 sec -1 , U c /kT = 12.62 and u = 2.03. 

Fig. 5 a) The field profiles B{x) for flux creep at U = U (B/ H*)(j c /j — 1) with H* — H/2 found by numerical solution 
of Eq. ( ]i"o| ) (squares) and by the semi- analytical approach (solid lines), 
b) U(x) found for the same U(B, j)-dependence from the numerical solution of Eq. (fio|). 

Fig. 6 a) The same as in Fig. 5a for the collective creep dependence U — Uq(B / H*) a (j c / j) M with a = 1, a = 2 and 
H* = H/2. 

b) Magnetization for the same U(B,j )-dependence found from the numerical solution of Eq. (circles) and 
by semi-analytical approach (solid line). For m/H < 0.1 the circles and the line completely coincide. 

Fig. 7 Dynamic development of anomalous magnetization (fishtail) found by the semi-analytical solution of the Kim 
model (see Eq. (|4C|)) with B c = 2H* , a = 1 and a = 2. Relaxation starts at t = from j c (H) = joB c / (B c + H) 
shown as dashed line. Due to faster relaxation at small H an anomalous magnetization develops at j <C 
j c . Circles and solid lines correspond to the direct numerical and semi-analytical solutions, respectively, for 
\n(t/T f i ow ) = 2.8, 7.4 and 14.3. 

Fig. 8 The magnetic induction B (squares), vortex velocity v (circles) and the magnetic current D — Bv (triangles) 
vs. d— x found numerically from Eq. ([ll]) near the sample edge x — d. The fits are: \J d — x for B and l/y/d — x 
for v. Note that D shows no peculiarity at x = d. 

Fig. 9 The activation energy U/kT in the remanent state: a) a — 0.1, fj, = 1; b) a = 1.5, n = 0.5 c) a = 1, fi = 1. 
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